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J> ' Abstract 

A Borg-Marchenko type uniqueness theorem (in terms of the Weyl 
function) is obtained here for the system auxiliary to the iV-wave 
qq . equation. A procedure to solve inverse problem is used for this pur- 

pose. The asymptotic condition on the Weyl function, under which 
the inverse problem is uniquely solvable, is completed by the new and 
simple sufficient condition on the potential, granting the fulfillment 
of this asymptotic condition. The evolution of the Weyl function is 
5_i ■ discussed and the solution of an initial-boundary value problem for 



the iV-wave equation follows. Explicit solutions of the system are 
obtained. System with a shifted argument is treated. 

1 Introduction 

The well-known integrable nonlinear optics (iV-wave) equation has the form 
[D, u t \ - [D, u x ] = \[D, u], [D, u]\ , [D, u t ] := Du t - u t D, (1.1) 
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where D = D* and D = D* are m x m diagonal matrices, u x and u t are the 
partial derivatives of the m x m matrix function u(x, t) = u(x,t)*. Equation 
( 11. ID is the compatibility condition of the auxiliary systems 

Y x = G(x,t,z)Y, Y t = G(x,t,z)Y, 

where 

G:=izD-C(x,t), G:=izD-((x,t), £=[D,u], (=[D,u\. (1.2) 

The case of three waves (n = 3) interaction have been treated in the seminal 
paper [U], using zero curvature representation. It was one of the first models 
to demonstrate the advantages of the zero curvature representation. (See 
also [3] for the iV-wave case.) Several years later model (11.11) proved to be 
the first integrable system, which was studied via Riemann-Hilbert problem 
approach [JT]. The iV-wave equation is actively studied and used in optics, 
fluid dynamics and plasma physics. It describes a variety of phenomena, 
including patterns and various instabilities, and is also closely related to 
many other important nonlinear integrable wave equations. For these and 
other connections and applications see, for instance, the books and recent 
publications [IlEllSElESlESESlEIlEEland references therein. 
We shall consider the auxiliary to (II. ip system on the semiaxis: 

Y x (x,z) = (izD-£(x))Y(x,z), x>0 (Y x :=-^), (1.3) 
where, without loss of generality, we assume D > 0: 

D = diag {d l} d 2 , ■ ■ ■ , d m }, d x > d 2 > . . . > d m > 0; ((x) = -C*(x). 

(1.4) 

Here diag stands for diagonal matrix, and ((x) is an m x m matrix function. 
The study of system (II. 3p is basic for study and construction of the solutions 
of the iV-wave equation (for study of the difficult and important initial- 
boundary value problem for the iV-wave equation, in particular). Moreover, 
this system is a natural generalization of the well-known Dirac (called also 
Zakharov-Shabat or AKNS) system. Various useful results and references on 
the scattering problem for system (11.31) one can find in [21 El El El El ESI EZ] ■ 
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An extensive amount of research on the scattering problem for system (11.31) 
with the complex-valued entries dk of D was done by R. Beals and R.R. 
Coifman. The unique solvability of the problem on a suitable dense set of 
scattering data was obtained, in particular. The case (11.41) of the positive D, 
that is considered here, is less general, and therefore more explicit description 
of the class of data for which the solution of the inverse problem exists and 
a close to the classical procedure for solving this inverse problem have been 
obtained in [301 ED E2] in terms of the generalized Weyl functions. Notice 
that various generalizations of the Weyl functions are successfully used both 
in the inverse scattering and in the inverse spectral problems. 

Here we develop further the results from [301 ED EH] • It is shown that the 
asymptotic conditions (I2.9p and (12.101) on the generalized Weyl function, un- 
der which the inverse problem is uniquely solvable, are automatically fulfilled 
for the integrable and two times differentiable ( with integrable derivatives. 
If conditions (12. 9p and (12.101) are valid for the Weyl function of the initial 
system, they are valid also for systems with a shifted argument, and they 
hold under the evolution of the Weyl function too. In this way we obtain 
sufficient conditions, when the inverse problem for system (II. 3p and an initial- 
boundary value problem in the quarterplane for the iV-wave equation (11.11) 
are solvable, and our procedure for solving these problems works. 

A recent series of papers by F. Gesztesy, B. Simon and coauthors on the 
high energy asymptotics of the Weyl functions and local Borg-Marchenko 
type uniqueness results has initiated a growing interest interest in this impor- 
tant domain (see [121 [201 122l EHl H21 H3] and references therein). The Weyl- 
Titchmarsh theory for a non-self-adjoint case (the skew-self-adjoint Dirac 
type system) has been studied in [IT], [231 EH] and the Borg-Marchenko type 
results for this system have been published in [36]. In this paper we obtain 
a Borg-Marchenko type theorem for another important non-self-adjoint case, 
that is, system (11.31) . 

Construction of the explicit solutions of the inverse problems and non- 
linear equations is of great interest and Backlund-Darboux transformation is 
one of the most fruitful methods to do it. The initial approach by Backlund 
and Darboux have been greatly generalized and developed (see, for instance, 
PH ESI EH1 ED ESI US] and references therein). We consider some applications 
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and developments of the version of the Backlund-Darboux transformation, 
which is called GBDT (see [321 El [37] ) . 

Section 1 is an introduction and Section 2 contains preliminaries to make 
the paper self-sufficient. Section 3 contains Theorem 13.11 which states that 
for the two times differentiable ( with integrable derivatives asymptotic con- 
ditions (12.91) and (I2.10p on the Weyl functions are fulfilled. It contains also 
a Borg-Marchenko type Theorem 13.21 Weyl functions for systems with a 
shifted argument are treated in Section 4. Evolution of the Weyl function 
and solution of the initial-boundary value problem are described in Remark 
14.91 GBDT for system (II. 3p is discussed in Section 5. Lemma I4T71 is proved 
in Appendix A and formula (15.61) is proved in Appendix B. 

The space L 1 (a, b) = L^ xm (a, b) of m x m matrix functions on (a, b) 

is equipped with the norm ||/||i = J b \\f(x)\\dx, where the matrix norm is 

i 

defined in terms of the trace Tr by the equality ||/(x)|| = (Tt (f (x)* f (x))^j 2 . 
We denote by C k (a, b) the class of k times differentiable matrix functions. 

2 Preliminaries 

We shall consider system (11.31) . such that the inequalities 

sup ||C(^)|| < oo (2.1) 

0<x<l 

are true for each I < oo. The m x m fundamental solution w of system (11.31) 
is normalized by the condition 

w{0,z)=I m , (2.2) 

where I m is the m x m identity matrix. A generalized Weyl function, later 
called a Weyl function, is introduced for (II. 3p slightly different from [33] . 

Definition 2.1 A Weyl function of system U.3\) is anmxm matrix function 
ip, such that for some M > it is analytic in a lower semiplane < —M, 
and the inequalities 

sup z)(p(z) exp (— ixzD) || < oo (2.3) 

x < l, 3z< —M 

hold for all I < oo. 
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System (11.31) with a bounded on the semiaxis potential £: 

sup ||C(x)|| <M , (2.4) 

0<x<oo 

was treated in Lemma 1.1 [33]. For this case, it was proved that a Weyl 
function always exists and admits normalization 

<Pkj( z ) = 1 for A; = j, <p kj (z) = for k > j. (2.5) 

Moreover, by Theorems 1.1 and 2.1 [33J a normalized, as in (I2.5p . Weyl 
function of system (11.31) with a bounded on the semiaxis potential is unique. 
This Weyl function satisfies for some r > the inequality 

oo 

J ^exp (ixz D) S j(p(z)*w(x, z)*w(x, z)ip{z) exp ^x(— izD —rl m )j dx < oo, 
o 

(2.6) 

where Qz < —M. Inequality (12.61) is somewhat similar to the inequalities 
characteristic for the classical Weyl functions. 

Definition 2.2 The inverse spectral problem (ISpP) for system ( fl.gj) . which 
satisfies \2. is the problem to recover the system, i.e., to recover the matrix 
function (, such that 

C(x) = -C(x)*, ( kk (x) = 0, (2.7) 

from a Weyl function. We shall denote by Q the operator mapping the pair 
D and (f(z) into (, that is, Q(D,tp) = 

We no longer assume that <p(z) satisfies (12.51) . 

Theorem 2.3 [33] For any matrix function (p(z), which is analytic and 
bounded in the semiplane Qz < —M, and which has the property 

oo 

J {cp{z) - I m )* (cp{z) - I m ) d\<oo (z — X- ir), A g R, r } >M), 

— oo 

(2.8) 

there is at most one solution of the ISpP, i.e., Q(D,tp) is unique. 
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The existence of the ISpP solutions was proved in Theorem 1.3 [33] under 
stricter conditions. Namely, we require that 



sup 



z[(p(z) 



< oo (Ssz < -M), 



(2.9) 



and that for some matrix a for all lines z = X — ii] with fixed values rj > M, 
we have 



z^f{z)-I m -a/z) G L 2 mXm ( 
Without loss of generality we suppose also that 

det (p(z) ^ 0. 



-oo, OO . 



(2.10) 



(2.11) 



Theorem 2.4 [33] Let the analytic matrix junction ip satisfy A2.y\) - I[2~TT\) . 
Then ip is a Weyl function of a unique system M.3\) , such that \2. 7| ) holds. 

Finally, let us give the procedure to recover ( from ip (see the proof of The- 
orem 1.3 [33]). First notice that relations fl2.9l) - fl2.11l) yield 



sup 



z^z)- 1 -I m ^ <oo {Qz<-M), 



and 



a 



/z)eL 



2 ( 

mxin\ 



OO, OO ) . 



(2.12) 
(2.13) 



Therefore, we can introduce the m x m matrix function 



oo 

i r i 

IL(x) = / -(exp(ixzD)) (piz)' 1 dX 

2wi J z 



(z = X—ir), t}> M, x > 0), 

(2.14) 

where the integral is understood as the matrix function, which entries are 
the norm limits in L 2 (0, 1) of the integrals from a to b (a — ► — oo, b — > oo) of 
the entries of z -1 ( exp {ixzD)^j p(z)^ 1 . So H(x) G L^ xm (0,Z) is defined on 
each interval (0,/). Recall that 

oo 

/ exp (ixzD)/z dX = I m (x > 0). (2.15) 

2m J 
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According to (I2.12j) - (j2.15j) . U(x) is twice differentiable and the following 
properties hold: 

n(0) = J m , IL'(0) = -iDa, e- xMD W(x)eL 2 mxm (0,oo) (W = U x ), 

(2.16) 

e- lMD n"(x)eLi xra (0,oo). (2.17) 
Now, let us introduce the linear operator Si, which is bounded on L^(0,Z): 

i 

S l f = D- 1 f + J s(x,u)f(u)du, (2.18) 
o 

where s(x,u) = {s k j(x, u)}™ ; < x, u < I; 

7 

, c?r 1 IIl.„-(a; — djdl 1 !!,) for u<d k d~ 1 x, 
+ <( fc J fc (2.19) 

dJ x Hj k {u — d k d~ x x) for d k d~ x x < u; 

9(x,u) = {e kj (x,u)}l j=1 = n>(x)[n'(u)}*D-\ 

7 is the interval [max (0, x — d^d k u), x]. Sometimes we omit "/" in S l and 
write just S. The operator S satisfies the operator identity 

AS-SA* = ill IT, (2.20) 

X 

where the operator A acts in L^(0,Z): Af = iD J f(u)du, and II acts from 

o 

C m into 1/^(0, 1): Ug = U(x)g. Using the operator identity, it was shown in 
[33] that 

S t > e(l)I (e > 0, I - identity operator). (2.21) 
By (I2.16l) - fl2.21l) one gets the representation [2] 



X 

S-i = y*V, Vf = D^f + J V(x,u)f(u)du, (2.22) 



o 

7 



I X 

where J J V(x, u)*V(x, u)dudx < oo. The operator V' 1 admits representa- 
o o 

tion 

X 

V~ l f = D^f + J u)f(u)du. (2.23) 



Moreover, according to formula (1.60) in [33] we get 

V(u,x)* = -S- 1 s(x,u)D 1 ' 2 (0<x<u), 

u 

F(x,u) = s(x,u)D 1/2 + s(x,v) V(u,v)*dv (x>u>0), 



[2.24) 



o 



where S u 1 is applied to s(x,u) columnwise. Formula (I2.24p implies 

T(l, I) = D- 1 (Sfhfa I)) (I) D 1 ' 2 . (2.25) 

The potential ( can be easily recovered from r(Z, I) [33] : 

C(i) = (r(i,i) -DYil^D-^D 1 ' 2 . (2.26) 

Thus, formulas flgZED , fl^T8|) . fT2TT9|) . fl2T25|) . and fl2T26|) determine C The 
proof of (12.261) in [33J is based on the representation of the fundamental 
solution w: 

w(x, z) = P(x)w A (x, z), (2.27) 
where w A is the transfer matrix function in the Lev Sakhnovich form [55] - [40] : 

w A (l, z) = I m + izITS- 1 (I - zA) _1 n, (2.28) 

and f3(x) = (VU)(x). By the last relation, in view of (12722]) and (127241) . we 
have 

(3(1) = £>i(n(l) - (nOO.S-yt*,/))), (2.29) 

where (•, •) denotes a matrix with the entries, which are scalar products of 
the columns of matrices in the parenthesis, that is 

(n^S-yu,/))) =J (S , - 1 s(tt,Z))*n(«)dtt. (2.30) 
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Now, consider system (II. 3p on the whole axis, and let m x m matrix func- 
tion W satisfy (II. 3p . Then, matrix function A4(x, z) = W(x, z) exp(-ixzD) 
satisfies equation 



M.Jx, z) = iz 



D,M(x,z) -((x)M(x,z), -oo < x < oo, (2.31) 



and vice versa. The function Ai is defined by (12.311) up to the right factor 
exp(ixzD)a(z) exp(— ixzD). Normalization conditions 

lim M(x,z) = I m , lim a ._ >00 ||.M(2;, z)\\ < oo (2.32) 

are used to define M. uniquely. Sufficient conditions on the potential (, under 
which relations (12. 9p and (I2.10p hold, follow from a particular case of the very 
useful Theorem 6.1 [6j: 

Theorem 2.5 Suppose that themxm potential ( is two times differentiable, 
i.e. C(x) G C 2 (-oo, oo), and that (( k) {x) = e L^-oo, oo) for k = 
0,1,2. Then, for some M > the analytic in z matrix function Ai(x,z), 
which satisfies \2.31\) and 112. 32\) . is well defined in the domain Qz < —M, 
the norm is uniformly bounded: 

sup \\M{x,z)\\ < oo, (2.33) 

oo,oo), Srz<— M 

and uniformly with respect to x we have 

lim M(x,z)=I m . (2.34) 

\z\-^oo,$Hz<-M 

Finally, there is an m x m matrix function M.\(x) G C 2 (— 00, 00), such that 
\\M{x,z)-l m -—^-\\=0{z- 2 ) (9fc<-M), (2.35) 
and that M[ (x) G L 1 (— 00, 00) for k = 1,2. 

Scheme of the proof. First, the case ||£||i < 1 is treated in Theorem 3.8 
[B]. The operator K z ^ in the space of bounded matrix functions f(x) is 
introduced by the formula 

POO 

K z> cf{.x)= \ exp(i(x-y)zD)(C(y)f(y)) + exp(-i(x-y)zD)dy 

J X 
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exp (i(x - y)zD) (((y)f(y))_exp ( - i(x - y)zD)dy, (2.36) 



where g + (g_) is the upper (lower) triangular part of the matrix g, g = 
g + + g_, and the main diagonal is included in g_. One can see that ||i^ z ,^|| < 
||C ||i and so the operator I — K z ^ is invertible. Moreover, it is proved that 
A4(x,z) = (I — K z ^)~ l I m satisfies (12.311) and (12.321) . and all the properties 
of Ai follow. The general case \\C\\i < 2 r (r > 0) follows by induction on 
r (Theorem A |6J). In Theorem 6.1 matrix functions Aij(x) (j = 1,2) are 
introduced by the equalities lim^-oo M.j{x) = 0, 

-^Mj-i + CMj-x = i[D, M°% M%x) = -£ {ay)Mf(y)) d dy, 

(2.37) 

where Mo = I m , Mj = Mj + M° d , and the entries of Mj and Mj coincide 
on the main diagonal, the entries of Mj and M° d coincide outside the main 
diagonal. Then, relation (I2.35P follows from equality (6.14) in [6]: 

/"OO 

M (x, z)~ 1 M a (x, z) = I m — / exp(i(x — y)zD)g(y) + exp( — i(x — y)zD)dy 

J X 

+ / exp (i(x - y)zD)g(y)_exp (- i(x - y)zD)dy, (2.38) 

J — oo 

where M a (x, z) = Y^j=o z~*Mj(x) and 

g(x, z) = z~ 2 M'\x, z) (—M 2 (x) + C(x)M 2 (x)y 

3 Solvability of the inverse problem and Borg- 
Marchenko type result 

Theorem 12.51 yields our next theorem. 

Theorem 3.1 Suppose that the mxm potential ( is two times differentiable, 
i.e. CO) e C 2 [(D, oo), and that C (k) ( x ) = e ^t ; °°) f or k = 0,1,2. 
Then, for some M > a Weyl function of system $1.3\) satisfies conditions 

/ TO) and nnw . 
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Proof. Define ((x) on the semi-axis x < so that the conditions of Theorem 
12.51 hold. Then we have 

w(x, z) = M{x,z)exp{ixzD)M{0,z)- 1 (x > 0). (3.1) 

Hence, in view of Definition 12.11 and formulas (I2.33P and (13.11) . the function 
(p(z) = Ji4(0,z) is a Weyl function. Now, it is immediate from (I2.35j) that 
conditions fl23]) and fl2TT0D are fulfilled. ■ 

From the procedure to solve ISpP, a Borg-Marchenko type result follows. 

Theorem 3.2 Let the analytic m x m matrix functions tpi and if 2 satisfy 
i2.9i) - WHS) . Suppose that on some ray c^sz = $lz < f c 6 R, 3z < — M) 
we have 

ifiiz)- 1 - if 2 {z)' 1 = e~ ilzD 0{z) for \z\ -> 00. (3.2) 

Then ipx and if 2 ore Weyl functions of systems U.3\) with potentials (1 and 
(2, respectively, which satisfy {2. 7| ) and the additional equality 

Ci(x) = ( 2 (x) (0 < x < I). (3.3) 

Proof. The fact that ip\ and if 2 are Weyl functions follows from Theorem 
E2Q From m- flZHD follow relations f|27T2|) and fl2~T3|) . According to the 
classical results on the Fourier transform in the complex domain (see, for in- 
stance, Theorem V in [28]), the function -ip^z) -1 , where <f(z) satisfies (12.121) 
and (I2.13p . admits Fourier representation. Moreover, taking into account 
also formula (12.141) and Plansherel's theorem, we get this representation, for 
z = A — if) and fixed 77 > M, in terms of II: 

1 f°° 

-if(z)- 1 = iD (exp(-ixzD))U(x)dx, (3.4) 
z Jo 

where (exp(-xMD))ll(x) G L^ xm (0, 00). As we have ( exp(-xMD))ll(x) G 
L^ xm (0, 00), so equalities ( 13.41) hold pointwise. Hence, we can use (13 .41) to 
apply Phragmen-Lindelof theorem. Namely, put 

F(z) := (exp(ilzD)) (exp(-ixzD)) (ll^x) - U 2 (x)^jdx, (3.5) 
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where IIx and U 2 correspond via formula (12.141) to if\ and ip 2 , respectively. 
By (JS3D and (JX5D we obtain 



F(z) 



-D- l (exp(ilzD)) (ipxiz)- 1 - ^(z)- 1 



-J (exp(i(l-x)zD) S )(ll 1 (x) -U 2 (x)^dx } (3.6) 



In view of (13.21) . the relation 



-D- 1 (exp{ilzD)) Ui{z)- 1 - ^(z)- 1 



is true. Recall that z — A — irj, r\ > M, and that 
(exp(-xMD))ll k (x) e L^ xm (0, oo). It follows that 



0(1) 



(3.7) 



J (exp -x)zD)^j (ill 



x) — Il 2 (x) )dx 



T) -> oo. (3.8) 



According to f l 3 . 6 f> - ( T378T) . the matrix function F(z) is bounded on the ray 
c^z = IRz (Qz < —M). It is immediate that F is bounded also on the axis 
Qz = —M. Therefore, function F given by (13. 5p satisfies conditions of the 
Phragmen-Lindelof theorem in the angles with the boundaries Qz = —M 
and c^z = dtz (Qtz < —M) in the lower semiplane. That is, F is bounded 
for Qz < —M. It easily follows from (13.51) that F is bounded for Qz > —M 
too, and that F(z) — ► for z = z tending to infinity. So we derive F = 0. 
This implies that 

IIi(a;) =n 2 (a;) (0<x<Z). (3.9) 

Finally, notice that by (12.181) . (I2.19p . and (I2.25P the matrix function T(x, x) 
on the interval [0, I] is determined by U(x) on the same interval. Hence, 
formulas (l2T26|) and flM]) imply (Q. ■ 



4 System with a shifted argument 

In this section we shall consider system (II .3p with a shifted argument: 

Y x (x + <T,z) = (izD-((x + <r))Y(x + <r,z), x>0. (4.1) 
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Taking into account normalization condition (12.21) . for the fundamental so- 
lution w(x, a, z) of (14.11) we have 

w(x, a, z) = w(x + a, z)w(o~, z)~ . (4.2) 

By 04. 2 p and Definition 12.11 the next proposition is immediate. 

Proposition 4.1 Let f{z) be a Weyl function of system U.3\) . Then, the 
matrix function 

(p(a, z) = w(a, z)ip(z) exp(-iazD) (4.3) 
is a Weyl function of system 

Notice that according to (I2.24p the matrix function T(x,u) does not depend 
on the choice of / (/ > x > u > 0) for the domain of operators S and A. 
Putting in (12.241) x = u, we get 

X 

T^^sMD^ + Js^VM'^ (4.4) 



By f l2~T9|) . fl2T2TD . the first relation in fl2T24j) and equality ([PL one can see 
that T(x,x) is continuous for x > 0. (In fact, T(x,x) is differentiable.) 
Therefore, according to fl 2 . 2 6 [) the matrix function ((x) is continuous too. 
We shall put 

T(0,0) := lim T(x,x), (4.5) 

x— >+0 

and similar to (12.261) assume 

C(o) = ( r(o, o) - d r(o, o)d- 1 )d 1 ^ 2 . (4.6) 

Now, we can express ((0) in terms of the matrix a, which is defined by 
if via representation (12.131) . For that purpose introduce an m x m matrix 
5? = {ttfcj}fcj=i via the entries a>kj of a: 

otkj ■= az kj for/c < j; a kj := -a jk fork > j. (4.7) 
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Proposition 4.2 Let the analytic m x m matrix function ip satisfy A2. 9\) - 
Then, for ( = Q(D, ip) we have 



C(0) =i(Da-aD). 
Proof. In view of (12.191) and (14. 4p one obtains 



(4.8) 



r(0,0)J fe = dfcMjnjyCO) for k <j; (T(0,0)J fe . = dJ 1 d]U jk (0)fork> j. 

(4.9) 

Recall that according to (12.161) we have II' (0) = —iDa, and so formulas (14.71) 
and (14.91) imply that T(0, 0) = —iaDz. Hence, formula (14.81) follows from 

ID. ■ 

Example 4.3 Consider the simplest example m = 2 and 



C(*0 



-q 
q 



const, 



(4.10) 



where const means a constant matrix. Calculate eigenvalues and eigenvectors 
of izD — C, to get 



izD-( = T(z)A(z)T(z)- 1 , A(z) 



Ai(z) 
X 2 (z) 



(4.11) 



T(z) 



(4.12) 



1 (izd 2 - X 2 {z))/q 

(Xi(z) -izdx)/q 1 
where A& are the roots of equation 

(izdi - X)(izd 2 - A) + \q\ 2 = 0, i.e. (4.13) 

Ai, 2 = ~((di + d 2 )z ± v / (rfi-rf 2 ) 2 ^ + 4|g| 2 ). (4.14) 
In particular, we have 

X k (z)-tzd k = (-l) k+1 2t\q\ 2 (V K - rf 2 ) 2 ^ 2 + 4|g| 2 + (di - d 2 )z) ~\ (4.15) 
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In view of \J^.ll ) we get the fundamental solution 

w(x, z) = T(z) exp{xA(z)}T (z)~ x . 



(4.16) 



From Iji4-15\ ) and (i4.16\) it follows that the matrix function (p(z) = T(z) sat- 
isfies 112.3]) , and so this ip proves a Weyl function of system II 1.3]) with ( of 
the form fi^Tffl . Moreover, from {4. 12]) and l[4-15\) it follows that 



<p(z) = T{z) = I 2 + 



(di - d 2 )z 

Therefore conditions ll2.9]) - l[2~Tl]) are fulfilled. 



q 
q 



(4.17) 



Remark 4.4 Notice that Example \4 .3] describes the simplest case, where con- 
ditions 112.9]) and 112.10]) are fulfilled, but conditions of Theorem \3.1\ are not 
valid. 

Using Proposition 14.21 we obtain a similar result for a system with a shifted 
argument. 

Proposition 4.5 Let the analytic mxm matrix function (p(z) satisfy $2.9]) - 
112.11]) . where 

a* = -a. (4.18) 

Then the matrix function ip{a,z) also admits representation $2.9\) - [2~Tl\) . 
where M{a) = M + e for an arbitrary fixed e > 0, and where the matrix 
a(o~) is such that 

a(cr)* = -a{a). (4.19) 
Moreover, for Q — Q(D, (p) we have 

((a) = i(Da{a)-a(a)D). (4.20) 

Remark 4.6 If the conditions of Theorem \3.1\ are fulfilled, we have 

<p(<r,z) = M(a,z), a(a) = M 1 (a). (4.21) 

According to the first relation in {2.31 ) and equality A4 = I m , we get 

[D,M 1 (*)} = -iC((T)=i((<Ty, M° d (o-y = -M° d (a). (4.22) 
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From the second equality in \2.31 ) and the first equality in \^.22 ) it follows 
that the k-th diagonal entry of M.\ has the form 

Mf(a)) =-i f Y.( d 3-d*Y 1 \CkM\ 2 dy = -{M d l {<r)) , (4.23) 



where M k = {j G N | < j < m, j ^ k}. Formulas j4-22\j and \J^.23\j 
imply M.\ = — M.\, and so, using {J^.21% we derive a* = —a, a(a)* = 
—a(a). Thus, under conditions of Theorem \3.1\ equality ff^JW is true and 
the statement of Proposition \4-5\ follows from |2J/. Still, the conditions of 
Proposition \4 . 5\ are weaker than conditions of Theorem \3.1\ (recall Example 

m 



To prove Proposition 14.51 we shall need some preparations. In view of (12.181) 
we can present S^ 1 (S = Si) in the form 



S- l f = Tf = Df + J T(x,u)f(u)du. 



From ST = I, according to fl2~T81 and it follows that 



(4.24) 



s(x,u)D + D~ l T(x,u) + J s(x,v)T(v,u)dv = 0. (4.25) 

o 

In particular, for the fixed values of u we shall assume 

T(x, u) = -S^six, u)D, (4.26) 

and we shall define also T(x,u) pointwise by the formula 

i 

T M = - D s M D + Dj(s^ x) )"s Miv D. (4.27) 



Introduce now an m x m matrix function K(x) by the formula: 

i 

K(x) := IT 1 (s^Uj (x) = U(x) + D' 1 J T(x, u)U(u)du. (4.28) 
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From the indentity fl2T2TIj) . it follows that S~ 1 A — A*S~ 1 = iS^UWS' 1 , i.e., 
i i 

I m + IT 1 J T(x, v)dv + J T(v, ujdvD" 1 = K(x)K(u)*. (4.29) 

11 X 

Using (I4.27P and (I4.28p . it is shown in the Appendix 1 that K(x) is continuous 
and differentiable. The following lemma is also proved in the Appendix 1. 

Lemma 4.7 Let the analytic m x m matrix function (p satisfy li2.9\) - l[KTl\) . 
and Iji4.18[ ). Then, the relations 

z(w{l, z)<p(z) exp(-ilzD) -I m - -K{1) (V'(Z)) V 1 ) G L 2 mxm (-oo, oo), 

(4.30) 

and 

sup z (w(l, z)ip(z) exp(— ilzD) — J m j < oo (4-31) 
are true for < —M — e (for any e > 0). 

Proof of Proposition 14.51 By formula fl4.3[) and Lemma 14.71 we see that 
<p(l,z) satisfies f[2T9>f[2TTil . and 

a(l) = iK{l) (K'il^D- 1 . (4.32) 

Taking into account (14.271) we get, that T(l,u) is continuous in u at u — I. 
Thus, putting x = I and differentiating both sides of f)4.29p with respect to 
u at u = I, we derive 

K{1) (K'{1)^ * = —D~ l T{l, l)D-\ (4.33) 

Hence, according to (14.321) and (14.331) . we have 

a(l) = —iD~ 1 T(l,l)D~ 1 . (4.34) 

From (EDI), (jZHSD , and fTCTSl) it follows that s(l,l) = s(l,l)*. Therefore, 
formula (I4.27[) implies T(l,l) = T(l,l)*, and so, in view of (14. 34j) . we have 
(I4.19p . According to (I4.19p . the equality a(a) = a(cr) is true, where 

a kj (a) := a kj (a) fork < j; a kj (a) := -a jk (a) for A; > j. 

Now, as ip(a, z) satisfies ( 12.9p -( !2TTTl) . formula ( 14.20f) follows from Proposition 
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Remark 4.8 Notice that under conditions of Theorem 2.4 formulas {2. 14 ), 
A2.18\) , 112.19]) , K2. 25\) . and l\2. 26}) define a solution ( of the inverse problem 
even without the requirement > dj for k > j . That is, if D > and 
dk 7^ dj for k ytz j, then (, which is recovered from (f by the mentioned 
above formulas, satisfies {2. 7| ) and defines such a system that ([Op holds. By 
Theorem \2.3\ the solution of the inverse problem is unique. Recall that we 
denote this solution ( by Q(D,if). 

Consider now nonlinear optics (iV-wave) equation (11.11) . where u = u* and 
D satisfies the second relation in (11.41) . 

Remark 4.9 We assume for convenience that the entries of u on the main 
diagonal are identical zeros, i.e., u^k = 0. 

Suppose first that the Weyl function ip is bounded and satisfies Ii2.8\) . 
Then Wl^ . for the case 

D = diag {d-L, d 2 ,..., d m }, d x > d 2 > . . . > d m > 0, (4.35) 

the initial condition 

[D, u(x, 0)] = tt(D, cp) (4.36) 

defines at most one continuously differentiable solution u of U.l\) on the 
semi-band x > 0, > t > 0. 

Consider now the more general case, where D satisfies the second relation 
in \l-4^ , D > and dk ^ dj for k ^ j , but the inequalities in \4 ■ 3ty do not 
necessarily hold. Define the initial condition via Weyl function if by formula 
lji4-36] ) and define the boundary condition via the same ip: 

[D,u(0,t)] = tt(D,ip). (4.37) 

Assume that the analytic mxm matrix function if satisfies li2.9\) - l[2~Tl\) . and 
(gTZSp. Then, we have JBTj: 

a) The evolution of the Weyl function is given by the formula 

ip(t, z) = R(t, z)ip(z) exp(-izDt), (4.38) 

where 

dR(t, z) 



dz 



(izD -Q{b,<f)^)R(t,z), R(0,z)=I m . (4.39) 
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Moreover, the matrix functions <p(t,z) also satisfy conditions li2.9\) - l[KTl\) . 
and for some matrices a(t). 



b) The matrix function u(x,t) = u(x,t)* is well-defined by the relation 

[D,u{x,t))=n(D,(p(t,z)) (4.40) 

and satisfies nonlinear optics equation U.l\) and initial- boundary value con- 
ditions U-36\) and j4.31\j . 



5 System with a shifted argument and Dar- 
boux matrices 

From formula (11.31) it follows that (p(o~, z) given by (14.31) satisfies equation 

^Lfl = iz(p<p{a, z) - <p(a, z)D) - z). (5.1) 

The next proposition is also true. 

Proposition 5.1 Letip(a,z) satisfy equation ( 15. and let(p(0,z) be a Weyl 
function of system M.3\) . Then, the matrix functions ip{a,z) with fixed a > 
are Weyl functions of systems OTTp. 



Proof. By (15.1ft one can see that ip(a, z)e lzaD satisfies ( 14. ip . i.e., 

if (a, z)e izaD (p(0, z)' 1 = w(a, z). 

Hence, <f(a, z) has the form (14. 3p . and our proposition follows from Proposi- 
tion ED ■ 

Notice that equation (15.11) coincides with the definition of the Darboux 
matrix, which transforms solution e lzxD of the auxiliary to nonlinear optics 
equation system with a trivial potential Co = into solution of system (II. 3p . 
Therefore, by constructing Darboux matrices we obtain examples of the Weyl 
functions. We propose below two schemes to construct Darboux matrices. 
The first scheme is a particular case of the so called GBDT (see [321 [3U [37] 
and references therein). Namely, we shall introduce Darboux matrix as a 



19 



transfer matrix function (in Lev Sakhnovich form) with additional depen- 
dence on the variable x: 

w A (x, z) = I m - iU(x)*S(x)-\A - zl m )- l ll{x). (5.2) 

Distinct from formula (12.281) . both Tl(x) and S(x) in the GBDT method are 
differentiable matrix functions, where II is determined by the linear differen- 
tial system. Correspondingly, the factor U(x)* above means multiplication 
by the matrix adjoint to II(a;). The second scheme is also constructed in the 
spirit of the GBDT approach. 

Scheme 1. This scheme is precisely GBDT for system ( II .3p (see [30], [32]). 
First, we fix an integer n > 0, two n x n parameter matrices A and S(0), 
and n x m matrix 11(0) such that 

AS(0) - S{0)A* = TI(0)n(0)*. (5.3) 

Introduce H(x) and S(x) for x > by the equations 

u x = -iAUD + nc, s x = norr, (5.4) 

and put 

C(x) = ((x) - (Dn(x)*5 , (ar)- 1 n(x) - U(x)*S(x)- 1 U(x)D) . (5.5) 

Then, in the points of invertibility of S(x), the transfer matrix function wa 
satisfies [32] the equation 

-^-wa(x, z) = G(x, z)wa(x, z) — wa(x, z)G(x, z), (5.6) 
ax 

where 

G(x,z) = izD-((x), G(x,z)=izD-((x). (5.7) 
By (I5.3p and (15.41) we have also 

AS(x) - S(x)A* = TI(x)n(x)*. (5.8) 

Further we assume that S(0) > 0. Then, according to (15.41) . we have S(x) > 
0, and so S(x) is invertible. To make the paper self-sufficient we give the 
proof of formula (15.61) in Appendix 2. 
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Scheme 2. Fix an interval [a, b] and an m x m weight matrix func- 
tion p{t) > 0, which is bounded on this interval. Thus, the space L 2 m {p) 
with the scalar product (/, g) = J g(t)*p(i)f(t)dt is generated. Let £(x) be 
continuous, and define operators A and S(x) in L^(p) by the formulas 

Af = tf(t), S(x)f = cf + z f w{x,t ) Mx,y) P (y)f(y)dy, (5.9) 

J a t-y 

respectively. Here we take the principal value of the integral in ( 15.91) . Notice 
that -^w(x, t)*w(x, y) = —i(t — y)w(x,t)*Dw(x,y), i.e., 



(t — y) 1 w(x,t)*w(x,y) = (t — y) l — w(u,t)*Dw(u,y)du. (5.10) 

Jo 

Hence, we can rewrite the expression for S from (15.91) in the form 

rb rb rx 

S(x)f = cf+i (t-y)~ 1 p(y)f(y)dy+ / w(u,t)*Dw(u,y)dup{y)f(y)dy. 

J a J a Jo 

(5.11) 

As p is bounded and the operator J a (t — y)~ 1 »dy in L 2 m is bounded, so the 
operator J (t — y)~ 1 p(y)»dy, and also operator S(x), is bounded in L^p). 
Next, introduce operator H(x), acting from C m into L^(p), and operator 
Il(x)*: 

U(x)f = w(x,t)*f, n{x)*g= [ w(x,y)p(y)g(y)dy. (5.12) 

J a 

By (15. lip and (I5.12p we have the second equality in (15. 4p . that is, S x = 
n/^n* > 0. Now, choose c so that 

cl + i [ (t-y)- 1 p(y).dy>0. (5.13) 

J a 

As S x > 0, by (I5TT3D we have S(x) > 0. In view of (O]) . definition (1531) 
of A and definition (15. 12ft . the first equality in ( 15.41) is true too. Moreover, 
according to (15.91) the identity (15.81) holds. By (15. 4p and (15.81) the equation 
(15.61) is satisfied, where G and G are given via (I5.5P and (15. 7p . - see Appendix 
IB! Using Definition 12.11 Proposition 14.11 and formula ( 15. 6p one easily gets 
the next proposition. 
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Proposition 5.2 Let the Darboux matrix wa be defined via i\5.2\) . using 
Scheme 1 or Scheme 2. Then the normalized fundamental solution of the 
transformed system 

Y x (x, z) = G{x, z)Y(x, z) = (izD - ((x))Y(x, z) (5.14) 

is given by the formula 

w(x,z) = wa(x, z)w(x, z)wa(0, z)' 1 , (5.15) 

where w is the fundamental solution of the initial system A1.3\) . Suppose 
also that ip(z) is a Weyl function of system ( f 1.3\) . Then the matrix function 
wa{0, z)ip(z) is a Weyl function of the system \5.14\ ) an d 



(p(a, z) = wa(o~, z)w(a, z)ip(z) exp(-iazD) 



(5.16) 



is a Weyl function of the system Y x (x + a, z) = G(x + a, z)Y(x + a, z) with 
a shifted argument. 

Next we shall consider 2 simple examples of Scheme 1, including the case of 
A non-diagonal. 



Example 5.3 Let 

C(x) = o, n = 2, n(o) = [/ a h ••• f m ] (f k ec 2 ). 

It follows from the first relation in fl5.^| ) that 

n(x) = exp(— ixdxA)fi exp(—ixd 2 A)f2 . . . exp(—ixd m A)f ri 

It is immediate that 



(5.17) 



(5.18) 



exp(— ixdA) = diag{exp(— ixda), exp(-ixda)} for A = diagja, a}, 

(5.19) 



and 



exp(-ixdA) = exp(-ixda) I I 2 — ixd 



1 




for A 



a 1 



a 



(5.20) 
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Finally, assuming a ^ a, formula ( f5. 8\) implies 

s = {skj}lj=i = i{a - a) _1 nn* 
for the case A5.19\) . For the case A5.20\) the same formula implies 

s 22 =j(a-a)" 1 (nn*) 22 , si 2 = s^= (a - a) _1 ^(nn*) 12 - s 22 ), 

s u = ( a -a)- 1 (i(im*) u + s 12 -s 2i y 

Substitute explicit formulas for II and S, which are given above, into A5.2\) 
and ( f5. 5\) to obtain explicit formulas for wa and (. 

Notice that wa(&, z) admits representation 

w A {a,z) =I m + % -U(ayS(a)- 1 U(a)+0^, z -> oo. (5.21) 

From Proposition 15.21 and formula (15.2ip follows corollary. 

Corollary 5.4 Let the conditions of Proposition \5.2\ be fulfilled, and let a 
Weyl function if of the initial system U.3\) satisfy li2.9\) - l[2~Tl\) with the cor- 
responding matrix a. Then, for some M > 0, a Weyl function tp(0, z) of the 
transformed system \5.1$ satisfies formulas \2. 5j) -( TOI]) . where matrix a is 
substituted by the matrix a: 

a = a + m(0)*S(0)- 1 U(0). (5.22) 

From Propositions 14.51 and 15.21 and from formula (15.211) we get. 

Corollary 5.5 Let the conditions of Proposition \5.2\ be fulfilled, and let a 
Weyl function of the initial system il.3\) satisfy formulas li2.9\) - l[2~Tl\) and 
H4.18I) with the corresponding matrix a. Then, for some M > 0, a Weyl 
function (p(a, z) of the transformed system with a shifted by a argument sat- 
isfies formulas $2.9\) - f2~Tl\) . UTWf} , where matrix a(a) is substituted by the 
matrix a(a): 

a(a) = a(a) + tU(a)*S(a)- 1 U(a) 1 ((a) = i{Da{a) - a{a)D). (5.23) 
Formulas (15331) yield (1531) . 
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A Appendix 

In this Appendix we shall obtain some properties of the matrix function K(x) 
defined in (I4.28p . and using these properties we shall prove Lemma [4.71 

Lemma A.l Let the analytic mxm matrix function (p satisfy Ii2.9\) - I[KTT\) . 
and h4-18\ ). Then, K(x) is twice differentiable and satisfies equalities: 

I m -(U'(u),K(u))=K(0)*, (A.l) 

(z) _1 n"(w), - (k'^Yd- 1 + (l m - (W^^Kiu^D-'U'iO) = 0. 

(A.2) 

Proof. According to (I2.9I) - (I2.11I) and (I2.14p . the matrix function U(x) is 
two times differentiable. By (14.271) the matrix function T(x,u) + Ds(x,u)D 
is continuous in both variables. Then, in view of formula (I4.28P K(x) is 
continuous. Consider now formula (14.291) . It follows that K(l)K(l)* = I m . 
Hence, for the case u = I formula (14.291) yields 



K(x) 



[l m + j T(v, ^dvD- 1 ^ K(l). (A.3) 



It follows from (1A.3[) that K(x) is differentiable and 

K'{x) = -T(x,l)D- 1 K{l). (A.4) 

Put now in (I4.29P x = I, multiply both sides from the right by IT (it), and 
integrate the obtained expressions with respect to u from to /. We have 

n(Z) — 11(0) + D~ x j [ T(l,v)dvW{u)du = K{1) [ K{u)*W{u)du. (A.5) 
Jo Ju Jo 
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Recall that by (12.161) we have 11(0) = I m and change also the order of inte- 
gration in (1A.5j) . Then we derive 

n(Z) - I m + D- 1 [ T(l,v){U{v)- I m )dv = K{1) f K(u)*W(u)du. (A.6) 
Jo Jo 

Using f)4.28p . rewrite flA.61) in the form 

K(l) - I m - D' 1 I T(l,v)dv = K{1) [ K( U yW(u)du. (A.7) 
Jo Jo 

According to (14.291) we get 

Im + D- 1 [ T(l,v)dv = K(l)K(0)*. (A.8) 
Jo 

Recalling that K(l)K{l)* = I m and using (1A.7I) and (1A.8I) . we finally obtain 
f lA.ip . By fT4~T8D and the second relat ion in (I2.16P we have 

D^U'iO) = ^-^'(0))*. (A.9) 

Hence, according to (12.191) . s(x,u) is continuous, and so T(x,u) is continu- 
ous. Moreover, in view of (14.271) and (1A.4I) one can see that K is two times 
differentiable and the entries of K"(x) belong L 2 (0, I). Formula (12.191) yields 
also 

D _1 n'(x) = s(x,0). (A.10) 

Put again in (I4.29P x = I, multiply both sides from the right by D -1 !!"^), 
and integrate the obtained expressions with respect to u from to /: 

zr 1 (n'(Z) - n'(o)) + zr 1 f t(i, v )d- 1 (n» - u'(o))dv 

Jo 

= K(l)(^D' 1 U"(u),K(u)y 
So, taking into account (lA.lOp . we have 

D~ l (ds(1,0) + T(l,v)s(v,0)dv \ - (l m + T{l,v)dv\ D^U'^O) 
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= K(l){D- l Y^\u),K{u)Y (A.ll) 
By M . fl426|) . and fl429l) . we rewrite (EOT]) as 

— D~ x T(l, 0)D~ X — K(l)K*(0)D~ l H'(0) = K{1) (d^H"^), K{u)\ (A.12) 

Recall that S 1-1 = (S 1-1 )* and that T(x,u) is continuous. Hence, it follows 
that T(x, u) = T(u, x)* and, in particular, that T(Z, 0) = T(0, /)*. Therefore, 
using (IA.4I) and the equality K(l)K(l)* = I m , we obtain 

T(/,0) = -DK(1)(k'(0)Y. (A.13) 

From (TA~T2l and ([03]) we get 

^'(0))* J D- 1 -^*(0) J D- 1 n'(0) = (d^II"^), (A.14) 

Finally, flSZEp and (jATTijl imply ([Oil . ■ 



Proof of Lemma 14.71 For the proof of lemma we shall use representations 
(E2ZD and (12T28D . First consider expression (7 - zA)' 1 !! from (jZZHD - It is 
easy to see that 

((7 - zA) -1 Il) (x) = U(x) + izD J exp (i(x - u)zD)U(u)du. (A.15) 

From (13.41) and (1A.15I) . using integration by parts, we obtain 

((/ - zA)' 1 !^ (x) 
= Il(x) + (expixzD) (^p(z)^ 1 — izD exp ( — iuzD)Yl{u)di 



= (expixzD)ip(z)- 1 + -D^iU'ix) + / exp (i(x - u)zD)U"(u)di 

z \ J x 

(A.16) 

By f[2T28l) and flATT6l) we have 

^^(z) exp ( — ilzD) = <f(z) exp ( — ilzD) + izWS^ 1 ^ exp (z(x — l)zD) 

+ % -D- l (l\'(x) + J exp(i(x-u)zD)Tl"(u)du^ip(z)exp(-ilzD)y (A. 17) 
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To consider the asymptotics of the right-hand side of (IA.170 we take into 
account that Tl*S~ 1 D~ 1 acts as the operator K(x)*»dx. Then, using inte- 
gration by parts, we get 

izWS- 1 exp (i(x - l)zD) = K(l)* - K(0)* exp ( - ilzD) 

- [ K'(x)*exp(ixzD)dxexp(-ilzD). (A.18) 
Jo 

Use integration by parts again to rewrite ( 1A.18I) in the form 

iz^S' 1 exp (i(x - l)zD) = K(l)* + -K'iiyD- 1 - K(0)* exp ( - ilzD) 



i 

z 

where 



(k'($)*D- x exp ( - ilzD) + q(l, *)) , (A.19) 



q(l,z):= / K"(x)*D 1 exp (ixzD^dx exp ( — ilzD), 
Jo 

and so, for any e > we have 

sup \\q{l,z)\\ < oo; (A.20) 

3K-A/-E 

and for the lines with the fixed values of S z we have 

q(l,z) G L^ xm (-oo, oo) (Qz<-M-e, -oo < 3? z < oo). (A.21) 

Consider now two other terms on the right-hand side of flA.17[) and take into 
account the second relation in (12.161) and flA.ip as well as the asymptotics of 
if, that is, formulas (I2.9P and (12.101) to obtain 

cp(z) exp ( - ilzD) - WS^D^Tl'^ipiz) exp ( - ilzD) 

= K{0)* (I m + -D^U^O)) exp ( - ilzD) +-qi(l,z), (A.22) 

where q± satisfies f)A.20|) and flA.211) . Finally, from integration by parts and 
asymptotics of (p it follows that 

POO 

-WS^D- 1 / exp (i(x - u)zD)U"(u)du(p{z) exp ( - ilzD) 

J X 
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= -( / K(xyD~ 1 Tl"(x)dx + q 2 (l,z^ = - ((ir^'fa), K{u)) + q 2 (l,z)} , 

*/ 

(A.23) 

where g 2 satisfies ( 1A.20I) and flA.211) . In view of flA.II) and flA.21) the sum of the 
right-hand sides of rfAlQ]) . (fA~22l and (TA~23l) equals X(0*+(z/z)^'(0* i:rl + 
q 3 (l,z)/z, where (73 satisfies ( 1A.20I) and ( 1A.21I) . In other words we have 



w A (l, z)<p(z) exp ( - ilzD) = K(l)* + -K^iyD' 1 + -q 3 {l, z). (A.24) 



Further notice that in view of (I4.26P and equality T(x, u) = T(u, x)* we have 
(S-hfal))* = ( — T(m, = ~D- l T(l,u). (A.25) 

According to (I2T29D , fl424j) and fTA~25|) it follows that 

f3{l) = D^D- l (S- l Tl){l). (A.26) 
Compare (14T28D and fTA~26l) to get 

/?(/) = D3lf(0. (A.27) 
By fl2T27l) , fLA24l . and fTA~26]) we obtain (OOI) and ■ 



B Appendix 



Proof of formula (15. 6p . 

From ([53D and ([EH]) it follows that 



* o— 1 



its 



iDTl*A*S- 1 - C^S' 1 - IT S^TIDTI* S 



* o— 1 



(B.l) 



By ([EHD we have = S~ X A - iS^IUTS- 1 , and so formula flRTj) can 

be rewritten as 



its -1 ) = (Dn*^^ - irs^ixD - oirs -1 + idtfs^a 

= -CITS -1 + {DIPS' 1 A, (B.2) 
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where ( is defined in ( 15.51) . Now, from from the definition (15.21) of w A and 
formulas H5.4j) and (IB. 21) it follows that 

= -CK - J ™) + zDirs- 1 ^ - zij-'u - (w A - I m )(-C) 

-U*S- X {A- zI m )' l AUD. (B.3) 

Substitute A = (A — zl m ) + zl m into the second and fourth terms on the 
right-hand side to rewrite (1B.3I) as 

^ = (izD -Q{w A - I m ) - {w A - I m )(izD - Q + DW S- l Il-W S- 1 T1D. 

(B.4) 

Formulas (15.51) and (1B.4|) imply 



(IZ 



dx 

and (15.61) is immediate. 



D - Q(w A - I m ) - (w A - / m )(zzD - C) + C - C, (B.5) 
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